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Abstract. Several interesting correlations among Gamma Ray Bursts (GRB) prompt and afterglow 
properties have been found in the recent years. Some of these correlations have been proposed also 
to standardize GRB energetics to use them as standard candles in constraining the expansion history 
of the universe up to z > 6. However, given the still unexplained nature of most of these correlations, 
only the less scattered correlations can be used for constraining the cosmological parameters. The 
updated fipeak — Ey correlation is presented. Caveats of alternative methods of standardizing GRB 
energetics are discussed. 
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INTRODUCTION 

Different cosmological probes seem to point toward a consensus scenario characterized 
by a universe which is experiencing an accelerated expansion. Standard candles (SNIa, 
e-g- [E] )> standard rulers (clusters, CMB and BAO, e.g. [2]) and the angular power 
spectrum of the CMB (e.g. [3]) are used to probe the Hubble expansion flow up to 
redshifts 1-2 (by means of SNIa or clusters) or up to the epoch of recombination (through 
the CMB anisotropics). 

Gamma Ray Bursts are observed up to very high redshifts: the farthest is 
GRB080913C at z=6.7 [4]. Among GRBs with known redshifts, 45% are at z > 2 
and 8% at z > 4. The high GRB luminosity (10^^ erg) and their detection in the y-ray 
band makes them attractive as a potential and complementary cosmological tool to 
constrain the cosmological models at z > 2. However, the problem is that GRBs are 
all but standard candles isj]: their isotropic equivalent energetics and luminosities span 
3-4 orders of magnitudes. Similarly to SNIa, it has been proposed to use correlations 
between vario us p roperties of the prompt emission [6^, |7j and also of the afterglow 
emission 



vario us p ro; 



11] to standardize GRB energetics. 



STANDARDIZING GRB ENERGETICS 

Isotropic energies (luminosities) can be computed for GRBs with measured redshifts 
and well constrained spectral properties. The spectrum gives the bolometric fluence 
F (peak flux P) and then E^.o = ^TtdtizfF /{I + z) = AndtizfP). These two 
quantities are strongly correlated with the rest frame peak energy fipeak of the vFy 
spectrum [12, Tsll . 



We have updated the sample of bursts with known z and spectral parameters to Jan 
2009 (the last being GRB 090102 at z=1.547). These are 97 GRBs. The fit of both corre- 
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FIGURE 1. Correlation between the rest frame peak energy and the isotropic equivalent energy (open 
symbols - 97 GRBs) or the collimation corrected energy (filled symbols - 29 GRBs). The solid (dashed) 
line is the best fit of the fipeak — Ey (fipeak — £^iso) correlation. The shaded regions show the l,2,3cT 
dispersion (modeled as a Gaussian) of the data points around the best fit lines. Red open symbols are 
the GRBs with a jet break time measurement reported in the literature. 



lations with a powerlaw gives £'peak°^£'iso'^'^^'^'^ '^^ and £'peak°^^iso*^''^''^'' '^'^- However, due 
to the large scatter of the data points is extremely large (492 and 612 for 95 degrees 
of freedom for the ^peak — ^^iso and ^peak — Aso correlations, respectively). The scatter 
of the data points is defined by their distance from the best fit line. By modeling the 
scatter distribution with a Gaussian we find a logarithmic dispersion of o =0.23 dex and 
a=0.28 dex for the fipeak ~ ^^iso and fipeak ^ Aso correlation, respectively. These scatters 
are much larger than the statistical errors associated with the observables (<7Epg^^) = 0. 1, 
(o^iso) = 0-06, (olj^o) = 0.07. There is then the possibility that a third variable is respon- 
sible for this large scatter. 

The large dispersion of the ^peak — £^iso and ^peak — ^iso correlations prevents their use 
to standardize GRB energetics. However, GRBs are thought to be coUimated sources. 
In the standard GRB model the jetted outflow should produce a break in the afterglow 
light curve decay [jl4|] . The measure of the time of occurrence of this break ?jet allows to 
infer the jet opening angle 0jet (under the standard afterglow model assumptions) and to 
recover the collimation corrected energy Ey = £'iso( 1 — cos 0jet) . 

By collecting available estimates of jet break times from the literature we found ^ 
and later confirmed jlSL \\6\ that the collimation corrected energy is strongly correlated 




FIGURE 2. Constraints on the cosmological parameters obtained with the fipeak — Ey correlatio n up - 
dated to Jan 2009 (29 bursts - solid line) compared to the previous update (19 GRBs - dashed line - flol '). 
Also shown are the constraints obtained from 156 SN type la (blue thin line - 117(1 ) and those from the 
WMAP data (green thin hne). 



with £peak- 

In FigH] the 

£^peak ~ Ey correlation is shown through the 29 bursts having 
a jet break in their optical afterglow light curves. We find fipeak = Ey'^^^^-'^^ with 

;t^=37.8 for 27 degrees of freedom assuming a wind-like profile of the circumburst 
medium. This correlation has a dispersion o = 0.07 which is consistent with the average 
statistical uncertainties on ^peak ^nd Ey and it is much smaller than the dispersion of the 
£'peak — £^iso and £'peak ~ Aso Correlations. This result also suggests that the dispersion of 
the ^peak — £^iso Correlation is due to the jet opening angle: by correcting Ei^o for Ojet (for 
each burst), the scatter of the fipeak — Eiso correlation is reduced. 

Due to its tightness the ^peak — Ey correlation can be used to constrain the cosmo- 
logical parameters (§, ^id]. The fipeak — Ey correlation, shown in FigUl was found by 
assuming a cosmological model (i.e. Q.m = 0.3, Q.^ = h = 0.7). It would be a circular 
argument to use this particular correlation to constrain the cosmological parameters. We 
originally solved this problem by properly accounting for the dependence of the corre- 
lation from the cosmological parameters dSJ or by adopting a Bayesian fitting method 

Figj2] shows the cosmological constraints obtained through the fipeak — Ey correlation 



with the most updated sample of 29 GRB. These are compared with the constraints 
obtained with the sample of 19 GRBs (l^. Also in this case we have applied the 
Bayesian method that overcomes the circularity problem. 



ISSUES 

The 

^peak ~ correlation is derived in the standard uniform jet scenario assuming a 
constant radiative efficiency and either a uniform [8] or a wind [15] circumburst density 
profile. UO] discovered a completely empirical correlation £[^0°^ ^peak^Jel between the 
three observables which are combined in the ^peak — Ey correlation. Also the empirical 
correlation, due to its low scatter, can be used to derive constraints on the cosmological 
parameters (l^. 

The mostly debated issue of the spectral-energy correlations in general is that they are 
due to selection effects flS]. It has also been claimed that the slope and normalization 
of the 

Epeak ~~ E'lso Correlation evolve with redshift lil9ll Selection effects can be studied 
in the observational plane corresponding to these correlations: the fipeak.obs — F and the 
£'peak,obs " P planes, where F and P are the bolometic fluence and peak flux. We have 
studied 21] two instrumental selection effects: the minimum flux required to trigger 



a burst (trigger threshold) and the minimum flux to properly analyze its prompt emission 
spectrum (spectral threshold). The former has been claimed to bias the £^peak~£^iso 
correlation [18]. The sample of 76 bursts (updated to Sept. 2007) with measured redshifts 
is composed by GRBs detected by different instruments. For this reason we modeled the 
trigger threshold and the spectral analysis threshold of the different detectors. 

We can exclude that the fipeak ~ ^^iso correlation is biased by the trigger threshold. 
Also we can exclude that the spectral threshold is biasing the pre-Swift sample. Instead, 
the Swift GRB sample (27 events) is biased by the spectral threshold. This is also due 
to the limited spectral energy range of the BAT instrument on-board Swift which limits 
the measure of ^peak.obs in the 15-150 keV energy range. By considering sub-samples 
of bursts at different redshifts, we also exclude that the ^peak ~ ^iso correlation slope or 
normalization change. 

Recently [1210 added GRBs without redshift in the £'peak,obs ~ F and £'peak,obs ~ P 
planes to compare the distribution of bursts with respect to the two selection effects. To 
this purpose a sample of 100 faint BATSE bursts, representative of a larger population of 
1000 objects, was analyzed. By means of this complete, fluence-limited, GRB sample, 
it was found that the fainter BATSE bursts have smaller fipeak.obs than those of bright 
events. As a consequence, the ^peak.obs of these bursts is correlated with the fluence, 
though with a slope flatter than that defined by bursts with z. Selection effects, which 
are present, are not responsible for the existence of such a correlation. About six per 
cent of these bursts are surely outliers of the fipeak — £^iso correlation, since they are 
inconsistent with it for any redshift. fipeak.obs also correlates with the peak flux, with a 
slope similar to the ^peak — Aso correlation. In this case, there is only one sure outlier. 
The scatter of the £'peak,obs — P correlation defined by the BATSE bursts of this sample 
is significantly smaller than the £'peak,obs — F correlation of the same bursts, while for 
the bursts with known redshift the fipeak ~ ^^iso correlation is tighter than the fipeak ~ Aso 
one. Once a very large number of bursts with ^peak.obs and redshift will be available, we 



expect that the ^peak — ^iso correlation will be similar to that currently found, whereas it 
is very likely that the £'peak ~^iso correlation will become flatter and with a larger scatter. 

One of the main drawback of the fipeak — Ey correlation for cosmological use is that 
it still has few points: the original sample of 15 events ^ has only doubled since 2004. 
With the launch of Swift [22] several jet breaks were expected to be measured, especially 
in the X-ray band. Jet breaks should be achromatic because they are produced by a 
geometric effect. All the jet breaks used to compute Ej in the pre-Swift era were, instead, 
obtained from the optical light curves. Swift revealed a complex X-ray afterglow light 
curve: the early afterglow is often characterized by a steep decay followed by a shallow 
phase lasting thousands of seconds [23] . A characteristic break time is that ending the X- 
ray shallow phase. [24] showed that this time is not a jet break. A possible interpretation 
is that the shallow phase is produced by a long lasting central engine activity [|25ll as 
also supported by the presence of strong precursors, post-cursors, and X-ray flares in a 
sizable fraction of bursts. Often the X-ray and the optical afterglow light curves do not 
track one another, suggesting that they are two different emission components. 

We selected a sample of 33 Gamma Ray Bursts (GRBs) detected by Swift, with known 
redshift and optical extinction at the host frame [26]. The de-absorbed and K-corrected 
X-ray and optical rest frame light curves are modelled as the sum of two components: 
emission from the forward shock due to the interaction of a fireball with the circum- 
burst medium and an additional component, treated in a completely phenomenological 
way. The latter can be identified, among other possibilities, as "late prompt" emission 
produced by a long lived central engine with mechanisms similar to those responsible for 
the production of the "standard" early prompt radiation. We find a good agreement with 
the data, despite of their complexity and diversity. Our approach allows us to interpret the 
behaviour of the optical and X-ray afterglows in a coherent way, by a relatively simple 
scenario. Within this context it is possible to explain why sometimes no jet break is 
observed; why, even if a jet break is observed, it is often chromatic; why the steepening 
after the jet break time is often shallower than predicted. 



CAN WE USE THE £:peak - ^iso OR £:peak - ^iso CORRELATIONS 

FOR COSMOLOGY? 

The use of the £peak — Ey correlation for cosmology requires to measure the redshift, 
the prompt emission £'peak,obs and ijet. The latter is the most critical observable: the jet 
break is typically observed at 1-2 days after the trigger and the light curve needs to be 
sampled at much later epochs in order to infer /jet when the afterglow can be very dim, 
also because Swift detects higher redshift bursts than before [27]. 

Therefore, it is interesting to explore if other correlations can be employed to stan- 



dardize GRB energetics. Recent attempts 11281 . 1291] tried to use the ^peak — £^iso and 



^peak ~ Aso Correlations. In this cases, however, one has to take into account that these 
correlations are affected by a dispersion which is much larger than the statistical uncer- 
tainty on the data points. 

The scatter of these correlations is due to three terms. One is the statistical uncertainty 
in the measurements of the parameters agtat- A second contribution to the scatter comes 
from systematic errors Csys which could also have a physical origin but is difficult to 
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FIGURE 3. Test of the use of the /Speak ~ £^iso correlation for constraining the cosmological parameters. 
The extra scatter terms of Eq.l, i.e. (7^,, vs Q.m are normalized to the average value of the statistical error 
associate to the variable they are assigned to. Top panel: the extra-scatter term is assigned to iipeak and 
the results similar to those of [28] are found. A minimum (though weak) is found for Q.m ^ 0.15. The 
solid curve is found with the likelihood function of Eq. 1 , the dotted and dashed curved are found with the 
symmetric likelihood function of |30]. The dot-dashed line is found with the fitting method. Bottom 
panel: the extra-scatter term is assigned to E\^q. This is also the most obvious assumption as it is E\^o to 
depend on the cosmological parameters. No minimum is found in this case of the extra-scatter term for 
any value of £Im- 



model. A third contribution to the scatter could finally be due to the cosmological model 
Ccosmo: in the real cosmology its contribution should be minimized. In the fipeak — £^iso 
and E'neak — Aso correlations the last two terms are dominating the scatter. 

Il28ll ( ll29il ) proposed the use of the ^peak — ^^iso (^^peak — Aso) correlation to constrain 
the cosmological parameters. What is appealing is the possibility to use a correlation de- 
fined by a large GRB sample (much larger than that defining the fipeak ~ correlation) 
and to use only on prompt emission observables (fipeak and ^isoOr Liso). 

However, in both correlations the non-statistical scatter a^yg + cr^osmo need to be 
modeled. These terms (combined) are treated as a free parameter, i.e. the extra-scatter 
Cext which is assumed to have a Gaussian distribution equal for all the data points. 

To test the possibility of using the ^peak ~ ^iso correlation for cosmology, [28] fit 
the fipeak — £'iso Correlation in different "cosmologies" (a flat Universe is assumed) and 



derive, in function of Q.m, the best fit values of the free parameters, i.e. the slope m and 
the normalization q of the correlation and the extra-scatter term Cext- The correlation is 
fitted with the likelihood function: 



1 



{ji-mxi-qY 



(1) 



where Xi.yi, Ox^i-, Oyj are the data points with their statistical errors and cr^, Oy are the 
projection of the extra-scatter CText along the coordinate axes. 



12811 assume Y=£'peak» X=£'isoand set the extra-scatter term = 0, i.e. they give the 
(free) extra-scatter only to ^peak- They find that the extra-scatter Oy shows a minimum 
corresponding Q.m ~ 0. 1 . Therefore, they apply the standard procedure to derive con- 
straints on the cosmological models through the ^peak ^ ^iso correlation. 

However, in the fipeak — ^^iso (^^peak — Aso) correlations it is E[^o (Aso) that depends 
on the cosmological parameters (through the luminosity distance Jl(z|^m, ^A:^o))- 
Therefore, the extra-scatter term should be assigned to X=Ei^o, i.e. Oy = 0. 

For this reason we repeated the same test on the ^peak — ^iso and ^peak — Aso correla- 
tions and we do not find any minimum of the extra- scatter Ox when it is assigned to Ei^o 
which is actually the variable that depends on the cosmological parameters. We verified 
our results also by (a) adopting the symmetric likelihood function of [30] (having the 
term 1 in the numerator of the first term of Eq.l); (b) fitting with the least square 
method; (c) inverting the order of the fitting variables (i.e. setting Y=£'isoand X=£'peak)- 
Similar results are also found for the fipeak — Aso correlation. Our results are shown in 
FigS 

In the Epeak — Ey correlation (Fig[T]) the scatter of the data points is already consistent 
with the statistical errors associated with ^peak and Ey. The only residual scatter is due to 
the cosmological model. This is why, without assumptions on the nature and "normality" 
of the unknown extra-scatter term, the E^^^ — Ey correlation is preferable to standardize 
the GRB energetics. The ^peak ~ Ey correlation also proves that most of the scatter of 
the fipeak ~ £^iso Correlation is due to the jet opening angle, which is different from burst 
to burst. 

This shows (as also recently demonstrated by fi^) that the 

^peak ~ ^iso and £peak ~ 

Liso correlations cannot be used straightforwardly to constrain the cosmological param- 
eters due to the unknown nature of the extra-scatter they are affected by. 
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